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ABSTRACT 

The recent discovery of an eclipsing hierarchical triple system with two low-mass stars in a close orbit 
(KOI-126) by Carter et al. (2011) appeared to reinforce the evidence that theoretical stellar evolution 
models are not able to reproduce the observational mass-radius relation for low-mass stars. We present 
a set of stellar models for the three stars in the KOI-126 system that show excellent agreement with the 
observed radii. This agreement appears to be due to the equation of state implemented by our code. 
A significant dispersion in the observed mass-radius relation for fully convective stars is demonstrated; 
indicative of the influence of physics currently not incorporated in standard stellar evolution models. 
We also predict apsidal motion constants for the two M-dwarf companions. These values should be 
observationally determined to within 1% by the end of the Kepler mission. 
Subject headings: Binaries: eclipsing — Stars: evolution — Stars: low-mass — Starspots 
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1. INTRODUCTION 

Double-lined eclipsing binaries (hereafter DEBs) are 
powerful tools for testing stellar evolution models. A 
wealth of information can be gleaned from observations 
of DEBs, including precise masses and radii of the com- 
ponent stars along with the apsidal m otion of the system , 
if the orbit is sufficiently eccentric (see lTorres et al.ll2010l 
for a review). To date, observations of DEB systems 
with at least one low- mass component (below 0.8 M ) 
have painted a grim picture for stellar evolution models. 
The radii predicted by models are systematically 5 - 10% 
smaller than those determined from observations, while 
the effectiv e temperatures derived from models are 5% 
cooler (e.g.JRibas et al.ll200l iMorales et al.ll2008. 2009; 



iTorres et al. 20101 andreferences therein). 

The aforementioned discrepancies have been at- 
tributed to the effects of large scale magnetic fields 
(|RibasH2006t IChabrier et Id] 120071) . as most of the low- 
mass DEBs discovered are close binary systems with 
short orbital periods (generally less than 3 days). Tidal 
interactions act to spin-up the stars as they progress to- 
wards complete tidal synchronization. Short rotation 
periods result from this process, amplifying the stellar 
magnetic field strength. Strong magnetic fields inhibit 
the efficiency of convective energy transport and increase 
the total surface coverage of starspots, effectively low- 
ering the temperature at the stellar surface. In order 
to co nserve flux, the stellar radius is forced to inflate 
(e.g-lGough fc Tavlerill966l : iMullan fc MacDonaldll200H 
IChabrier et al.1120071). 

Recently. iCarter et al.l (|2011l henceforth Cll) discov- 
ered a triply eclipsing hierarchical triple system, KOI- 
126, in the Kepler data set which contains two low- 
mass stars (KOI-126 B, C). The two low-mass stars 
are in a tight 1.77 day orbit itself orbiting a more 
massive primary star (KOI-126 A) on a fairly eccen- 
tric path about every 34 days. The authors were able 
to derive fundamental parameters for all three stars: 
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M A = 1-347 ± 0.032 M Q , R A = 2.0254 ± 0.0098 R e , 
M B = 0.2413 ± 0.0030 M Q , R B = 0.2543 ± 0.0014 R & , 
M c = 0.2127 ± 0.0026 M , and Rp = 0.2318 ± 0.0013 
Rq . Spectroscopy of the more massive primary indicated 
it has an effective temperature of 5875 ± 100 K with a 
super-solar metallicity ([Fe/H] = +0.15 ± 0.08). Lastly, 
the relatively large eccentricity of the systems allowed 
the authors to place a weak constraint on the apsidal 
motion constant for the two low-mass stars, determining 
that it is below 0.6 at the 95% confidence level. 

Cll derived an age for the system of 4 ± 1 Gyr 
after fitting KOI-126 A to theoretical Y 2 isochrones 
(jDemarque et al.l 12004). When Cll compared the ob- 
served mass- r adius relation to 1, 2, 4, and 5 Gyr 
iBaraffe et all (fl998l . BCAH98) isochrones for [Fe/H] < 
0.0, the models were seen to under-predict the radius of 
the two low-mass stars by 2-5%. However, it is impor- 
tant to note that no super-solar metallicity isochrones are 
available for the BCAH98 models (in part due to limita- 
tions of their equation of state). Cll suggested that the 
observed radius discrepancy could be the result of both 
the super-solar metallicity and possible magnetic activity 
of the system. 

Here, we present the results of theoretical stellar 
modeling of the KOI-126 compo nents using the Dart- 
mouth Stellar Evolution Program (Chabo yer et al.ll2001h 
IBiork fc Chaboverl 12001 IDotter et al l 120081 . DSEPL a 
descendant of the Yal e Rotating Stellar Evolution Code 
(|Guenther et al.|[l992l) . We derive model radii consistent 
with observations as well as theoretical apsidal motion 
constants for the two low-mass stars. In §2 we describe 
the DSEP models utilized for this study while in §3 we 
present our primary results. Finally, in §4, we discuss 
the implications of this study with regards to the only 
other well known low-mass eclipsing system, CM Dra. 

2. MODELS 

We constructed individual stellar evolution models 
for each KOI-126 component along with a series of 
theoretical isochrones using DSEF0. The physics in- 
corporated in the models have been described previ- 

3 http:/ /stellar. dartmouth.edu/~modcls/ 
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ouslv (iChabover et all l200lt iBiork fe Chaboverl [20061: 
iDotter et al.l 120071 12008ft . but we shall provide a brief 
summary. 

Our models include the effects of helium a nd heavy el- 
ement diffusion following the prescription of iThoul et al.1 
(1994), though for fully convective stars, diffusion physics 
are unimportant. The opacities utili zed by DSEP are 
the h igh-temperature OPAL opacities (llglesias fc Rogers 
11996ft and low-temperature opacities of iFerguson et al. 
(2005). Surface boundary conditions were defined us- 
ing the PHOENIX model atmospheres (jHauschildt et al] 
fl999albh and were attached to the interior model at 
T = T e ff by interpolating in model atmosphere tables. 
Attaching the atmospheres to T = T e ff makes the value 
of the convective mix ing length used in t he atmosphere 
code inconsequential (Baraffc ct al. 1997). 

Above 0.8 M Q , DSEP uses a general ideal gas 
equation of state with the Debye-Hiickel correction 
(|Chabover fc Kimlll995l ). In the low -mass regime, DSEP 
employs the FreeEOSQ in the EOS4 configuration, se- 
lected for its treatment of arbitrary heavy element abun- 
dances and its inclusion of the molecule. FreeEOS 
has also been shown t o be yalid f or modeling stars more 
massive than 0.1 M Q (|Irwinll2007l ). Con vective core over- 
shoot is included using the method of iDemarque et al.l 
(2004). Rotation was not considered. 

The only modification made to the underlying physics 
in DSEP is related to the partial inhibition of element 
diffusion. We have introduced t urbulent diffusion as 
described by iRichard et al.l (|2005ft . Turbulent diffusion 
modifies the atomic diffusion coefficient and acts to ex- 
tend the mixing region below the convection zone. The 
magnitude of the turbulent diffusion coefficient is tied to 
an adjustable reference temperature, To, and varies with 
density via: 

r 1 ~ 3 

D T = uD He (T ) (1) 
[p( T o). 

where uj characterizes the relative strength of turbulent 
diffusion and _D# e (To) and p(Tq) are the helium diffusion 
coefficient and density at the prescribed referenc e tem- 
perature, respectively. |Proffitt fc Michaudl (|1991[ ) moti- 
vate the p~ 3 dependence in order to reproduce the solar 
beryllium abundance, which appears to be unchanged 
over time. Thus, any non-standard mixing in the Sun 
must be localized to a narrow region below the solar con- 
vection zone. We select u> = 40 and leave it fixed, in 
concordance with IRichard et al.l (|2005l ). The reference 
temperature used primarily in our models is To — 10 6 , 
which was found to best re produce the obse rved abun- 
dance trends of NGC 6397 (|Korn et al.ll2007ft . 

A solar calibration model was generated to determine 
the appropriate initial mass fractions of helium (Yi n i t ) 
and met als (Zj n j t ), given th e solar heavy element compo- 
sition of lGrevesse fc Sau val (|1998j) . as well as to calibrate 
the convective mixi ng length, (cvmlt = t/Hp). At the 
solar age (4.57 Gvr: iBahcall et al.l l2005) we were able to 
reproduce the solar radius, solar luminosity, radius of the 
convective boundary, and (Z/X)q with olmlt = 1-938, 
Y tmt = 0.27491, and Z lmt = 0.01884. All of the models 
utilized in this study were calculated using the solar cal- 

4 By Alan Irwin: http://freeeos.sourceforge.net 
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Fig. 1. — Stellar evolution tracks used to constrain the age of 
KOI-126 A. The dark band signifies the radius constraints imposed 
by the observations. 

ibration as a reference and were assumed to be coeval. 
Super-solar metallicity models with [Fc/H] = +0.15 were 
generated with Y init = 0.28419 and Z mlt = 0.02469. 

3. RESULTS 

3.1. Stellar Age 

The age of the system was constrained by evolving a 
1.347 M Q model, with [Fe/H] = +0.15 and solar cali- 
brated ccmlt, and matching the model radius with the 
observed radius of KOI-126 A (Figure 1). The age we 
derive for the system is 4.1 ± 0.6 Gyr, consistent with 
the age derived by Cll. Our uncertainty in the age is 
dominated by the uncertainty in the observed mass and 
metallicity of the system, with the observational uncer- 
tainty of the radius being of negligible importance. The 
effect of changing the reference temperature for turbu- 
lent diffusion was also considered (including removing it 
entirely) and was found to play a negligible role in the 
age determination. 

3.2. Mass-Radius Relation 

The primary results of this Letter are demonstrated 
in Figure 2. DSEP accurately reproduces the ob- 
served radius for each of the low- mass stars at 4.1 
Gyr with [Fe/H] = +0.15. We find that for masses 
M 2 = 0.2410 M Q and M 3 = 0.2130 M Q the pre- 
dicted radii from DSEP are R 2 = 0.2544 R Q and 
R3 = 0.2312 Rq, indicating a relative error between the 
model and observed radii of less than 0.3%. At solar 
metallicity, the models predict radii approximately 1% 
smaller than those predicted by the super-solar metallic- 
ity models. The predicted radii are robust. Artificially 
reducing the mixing length (olmlt = 1-00) and fitting 
the atmosphere to a deeper point in the stellar envelope 
(t = 100) both produced radius changes under 0.5%. 

Solar metallicity isochrones from DSEP display radii 
approximately 1% larger than radii predicted by 
BCAH98. The difference is likely a consequence of the 
equation of state (EOS) utilized by each group. Whereas 
the EOS used by BCAH98 is calculated for a pure hydro- 
gen/helium plasma, FreeEOS calculates the EOS for an 
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CM Dra 
KOI-126 i 
DSEP [Fe/H] = 0.0 
DSEP [Fe/H] = +0.15 
BCAH98 [Fe/H] = 0.0 



0.22 0.23 
Mass (M - ) 

Fig. 2. — Mass-radius relationship as defined by CM Dra (red 
- points) and KOI-126 (black — points). Overlaid are 4.1 Gyr 
theoretical isochrones from DSEP with [Fe/H] = 0.0 (black - solid) 
and [Fe/H] = 0.15 (blue - dash) and from BCAH98 with [Fe/H] = 
0.0 (red - dash-dot). 

arbitrary metal abundance. Our models can be more reli- 
ably calculated above solar metalli city. To test the effect 
of the EOS, we ran DSEP with the ISaumon et al.l (|1995l 
SCVH) EOS as an attempt to mimic the BCAH98 model 
radii predictions. Our models using the SCVH EOS pro- 
duced radii within 0.5% of the BCAH98 models at the 
same mass and composition, illustrating the importance 
of the EOS. 

3.3. Relative Fluxes 

The effective temperature of a star is another mea- 
sure with which to compare our models. While Cll 
found the effective temperature of KOI-126 A to be 
5,875 ± 100 K from spectroscopy (which our models 
match), they were not able to determine the effective 
temperatures of the two M-dwarfs. However, the dy- 
namical/photometric model utilized by Cll yielded the 
flux for each M-dwarf relative to the primary. The 
model results were: f B /fA = (3.26 ± 0.24) x 10~ 4 and 
Self a = (2-24 ± 0.48) x 10- 4 . 

To compare our stellar models, we derived a synthetic 
color- T e f f transformation for the Kepler bandpass using 
the spectral response function provided in the Ke pler In- 
strument Handbook (jVan Cleve fc Ca ldwell 20091) . With 
model fluxes determined using the PHOENIX model at- 
mospheres, it was possible to generate synthetic colors 
for each of the model components of KOI-126. We found 
fB/fA = 4.77 x 10- 4 and fc/f A = 3.71 x 10" 4 , a 6a 
and 3(7 difference, respectively. We also derived relative 
fluxes utilizing the empirical relation given on the Kepler 
Guest Observer websitc0. The empirical relation relies on 
the SDSS g and r magnitudes to derive the approximate 
Kepler magnitude. SDSS magnitudes for our models 
were calculated using a synthetic color-T e ^ f transforma- 
tion (|Dotter et al.ll200"8l) . allowing for an estimate of the 
Kepler magnitudes. This semi-empirical transformation 
yielded: f B /f A = 3.07 x 10~ 4 and f c /f A = 2.36 x 10~ 4 , 
consistent with the relative fluxes derived by Cll. 

5 See http://keplergo.arc.nasa.gov/CalibrationZeropoint.shtml 



In an effort to decrease the relative flux for each M- 
dwarf in the purely theoretical transformation, we var- 
ied the parameters of the primary star within the given 
observational constraints and changed the amount of 
convective core overshoot (CCO). When considering the 
largest mass (1.379 M ), lowest metallicity (+0.07), and 
a relatively high amount of CCO, the relative fluxes were 
reduced to within 2a of the accepted values. As with 
the age determination, the tight observational constraint 
imposed on the radius of the primary made the range of 
relevant radii insignificant in deriving our results. 

3.4. Apsidal Motion Constant 

By the end of the nominal Kepler mission, Cll pre- 
dict they will know the apsidal motion constant of the 
two low-mass stars with a relative precision of nearly 1%. 
From our low-mass stellar models, we were able to pre- 
dict an apsidal motion constant for each star, quantifying 
the degree to which the mass within the star is centrally 
concentrated . To do this we followed the prescription of 
IKopall (|1978D and solved Radau's equation with j = 2: 

r^ + ^(* + l) + *fa-l) = 6 (2) 

where 772(7*) is the value of a particular solution to 
Radau's equation related to the stellar deviation from 
sphericity, p(r) is the density of the plasma, and (p) is 
the average density. We used a 4 th order Runge-Kutta 
integration method to solve Radau's equation to obtain 
a particular solution at the surface of each star. Having 
determined 772 (i£), we were able to determine the apsidal 
motion constant via: 



U = 



4 + 2772 (R) 



(3) 



As a check on our apsidal motion constant integrator, we 
generated polytropic models characterized by the poly- 
tropic constant n = 1.0 and 1.5 and compared the resu lts 
of our code with the results of iBrooker fc Ollel ([1955). 

The apsidal motion constants derived from the inte- 
rior structure of our models for KOI-126 B and C are 
k 2 = 0.1499 and k 2 = 0.1512, respectively. This suggests 
our models are slightly less centrally condensed than an 
n = 1.5 polytrope which is characterized by k 2 = 0.1433. 
In fact, the run of density and pressure for our models 
indicate our models are best described by a polytrope 
with n ~ 1.45. A set of theoretical apsidal motion con- 
stants for fully convective stars, generated using DSEP, 
is given in Table [T] and may be compared to future ob- 
servations. Finally, we note the effects of rotation on the 
de rived k 2 yalues a re negligible. Employing the formulae 
of lStothers! lfl974h . we find that rotation affects our k 2 
values at the 0.02% level for stars rotating with a period 
of 1.7 days. 

4. DISCUSSION 

Low-mass stars below the fully convective boundary 
are a wonderful tool to test basic physics. Their low mass 
affords theorists a stable, long lived (> 10 11 yr) labora- 
tory with which to test the physics incorporated in the 
models. These fully convective stars have relatively sim- 
ple structures, and uncertainties in the opacities, surface 
boundary conditions, and the treatment of convection 
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TABLE 1 

APSIDAL MOTION CONSTANTS (K2) FOR FULLY CONVECTIVE STARS WITH [Fe/H] = -0.50, 0.0, AND +0.15. VALUES ARE QUOTED AT 7 

DIFFERENT AGES (iN Gyr). 



Mass (M Q ) 


[Fe/H] 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


10.0 


0.20 


-0.50 


0.1523 


0.1523 


0.1522 


0.1522 


0.1522 


0.1522 


0.1521 




0.00 


0.1521 


0.1521 


0.1521 


0.1521 


0.1521 


0.1520 


0.1519 




+0.15 


0.1518 


0.1517 


0.1517 


0.1517 


0.1517 


0.1517 


0.1516 


0.25 


-0.50 


0.1499 


0.1498 


0.1498 


0.1498 


0.1498 


0.1497 


0.1496 




0.00 


0.1496 


0.1496 


0.1496 


0.1496 


0.1495 


0.1495 


0.1494 




+0.15 


0.1496 


0.1496 


0.1496 


0.1495 


0.1495 


0.1495 


0.1493 


0.30 


-0.50 


0.1482 


0.1482 


0.1481 


0.1481 


0.1481 


0.1480 


0.1479 




0.00 


0.1480 


0.1480 


0.1479 


0.1479 


0.1478 


0.1478 


0.1477 




+0.15 


0.1479 


0.1479 


0.1479 


0.1478 


0.1478 


0.1477 


0.1476 



have relativel y small effect s on the predicted properties 
of the models (|Dotterll2007l § 3.2). Modelers are relieved 
of having to specify a free parameter since the internal 
structure is insensitive to the prescribed mixing length. 
The age of these stars is also of little consequence, as 
the radii of low-mass stars is nearly constant over their 
main sequence lifetime. These effects imply a unique 
mass-radius relation for stars of a given composition be- 
low the fully convective boundary. The discovery of two 
new data points in this regime indicate the mass-radius 
relation is not unique and there is significant dispersion 
amongst the data. 

Discrepancies are still seen b etween our models a nd the 
components of CM Draconis ()Morales et al J 120091 here- 
after CM Dra), as shown in Figure 2. Comparing KOI- 
126 and CM Dra, both have low- mass components with 
strikingly similar masses, estimated ages, and orbital pe- 
riods (1.77 d and 1.27 d, respectively). Both systems are 
thought to be tidally synchronized, though not circular- 
ized. Their similar characteristics suggest the M-dwarfs 
in both systems should have similar radii. Not only do 
they not, but the opposite relationship of what is theo- 
retically ex pected is observed. CM Dra has a sub-solar 
metallicity ()Morales et al J 120091 and references therein) 
and should therefore possess smaller radii than KOI-126. 

It is possible that the KOI-126 system has been caught 
in a period of inactivity, similar to a solar minimum. CM 
Dra appears to have undergone such a period of quies- 
cence in 2000. Morales et al. noted that no corrections 
due to starspots were needed in the analysis of the light 
curve data from that year. If KOI-126 B and C are in 
a magnetically quiescent state, it is possible that their 
radii would not show signs of starspots or inflated radii. 
We note that starspots would likely be of negligible im- 
portance since variations due to spots would be reduced 
to noise amongst the signal of KOI-126 A. 

We find it encouraging that DSEP is able to predict 
the mass-radius relation suggested by KOI-126 B and 
C. Although, more data points are needed to allow for a 
more complete understanding of the reliance of the mass- 
radius relation on physics (standard and non-standard) 
incorporated in current models. It is clear from this work 
that no standard stellar evolution model will be able to 
simultaneously fit both CM Dra and KOI-126 and that 
work on non-standard stellar evolution models will be re- 
quired to fit CM Dra. As suggested by numerous authors, 
magnetic activity is likely the culprit and must be incor- 
porated into the next generation of models. However, 
the effects of a magnetic field on the interior structure 



of stars has previously been considered as a reduction in 
the prescribed mixing length, mimicking the reduction 
in convective efficiency that should accompany the pres- 
ence of a magnetic field. With this in mind, it is not 
clear how magnetic activity would affect stars below the 
fully convective boundary. Self-consistent magnetic stel- 
lar models should help lend insight into the discrepancies 
with CM Dra. 

Concerning the relative flux discrepancies observed be- 
tween the purely theoretical transformation and the pho- 
tometric models of Cll, the color- T e ff transformations 
are of interest and might not provide a fully accurate 
transformation to the observational plane. In the low- 
mass regime, opacities are complicated by the formation 
of molecules and the peak of the stellar spectrum is near 
the cut-off of the Kepler response function. Interest- 
ingly, the semi-empirical transformation yields relative 
fluxes which are entirely consistent with the photometric 
models of Cll. However, we must be cautious with this 
result as the systematic errors are not well constrained 
for the transformation from the SDSS magnitudes to the 
Kepler magnitude. It must also be noted that the flux of 
the primary star is sensitive to the details of CCO. We 
observed that increasing the amount of CCO brought the 
purely theoretical fluxes closer to the observational val- 
ues. More investigation will be required to accurately 
diagnose the discrepancies. 

Finally, we note that the determination of the apsidal 
constant will provide a crucial test of our stellar evolution 
models. In particular, it will test the EOS, which directly 
determines the run of density necessary for the computa- 
tion of the apsidal motion constant. Morales et al. found 
that using the BCAH98 models, was approximately 
0.11 whereas our models predict a larger value of approx- 
imately 0.15. The difference between these two values is 
directly attributable to the EOS, which determines the 
run of density within a stellar model. If Cll are able to 
accurately derive the apsidal motion constant to within 
1%, it will provide a stringent benchmark against which 
to test the interior physics of low-mass stellar evolution 
models. 

5. SUMMARY 

In their discovery paper, Cll reported that the triply 
eclipsing hierarchical triple KOI-126 appeared to sup- 
port the mounting evidence that current standard low- 
mass stellar models are unable to reproduce the observed 
mass-radius relation. However, we have generated stellar 
models and theoretical isochrone tracks using the Dart- 
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mouth Stellar Evolution Program and find our model 
radii agree with the observations. Combining the KOI- 
126 measurements with previous observations of the low- 
mass binary system CM Dra, we find that the dispersion 
in the observed fully convective mass-radius relation is 
significant and stands in contrast to theoretical predic- 
tions. The fact that CM Dra, a system with sub-solar 
metallicity, lies on the super-solar side of the theoreti- 
cal mass-radius relation is indicative of physics currently 
not incorporated in standard stellar models. We predict 
the apsidal motion constant for each of the KOI-126 low- 
mass stars and find &2 — 0.15. Cll postulate that they 
will be able to determine the apsidal motion constant 



with a relative precision of 1% by the end of the nominal 
Kepler mission. This will provide a crucial test for our 
models and will provide a stringent constraint against 
which to test all current and future low-mass standard 
stellar evolution models. 
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